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Spin Hall effect and Zitterbewegung in an electron waveguide 
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We study spin-resolved probability distributions for electrons in a multi-channel waveguide in the 
presence of spin-orbit interaction. For a spin-polarized electron injection, a Zitterbewegung pattern 
is predicted in the probability distribution of the electrons in the waveguide. For a spin-unpolarized 
injection, the spin-resolved electron probability in the waveguide shows spin accumulations. In 
addition to the spin Hall phenomenon, namely accumulations of opposite spins at the lateral edges 
of the waveguide, we predict the existence of a regular stripe pattern of spin accumulations in the 
internal region of the waveguide. We show that the predicted Zitterbewegung and spin Hall effect 
stem from the same mechanism and are formed from coherent states of electrons in the waveguide. 

PACS numbers: 72.25.Dc, 71.70.Ej, 85.75.Nn, 73.23.Ad 
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Spin-orbit interaction (SOI) is a relativistic effect: it 
arises from the fact that an electron moving with a fi- 
nite momentum in a non-uniform electric potential will 
see an effective magnetic field acting on its spin. The 
SOI plays an important role in the field of spintron- 
ics since it is in most cases responsible for spin rota- 
tion and relaxation. Recently an increased interest has 
been focused on a consequence of the SOI, namely a 
spatial spin accumulation in a two-dimensional system 
in the absence of a magnetic field, known as the spin 
Hall effect^^- 7 ! 8 ! 9 ! 10 ^ 1 ! 12 ! 13 ! 14 ! 15 ! 16 .!! Another SOI- 
induced phenomenon of current interest is the Zitterbe- 
wegung of electrons in semiconductors ) 18 ' 19 ! 20 in analogy 
with the oscillating behavior of a free relativistic electron 
due to interference between its positive- and negative- 
energy state components. The previous theoretical stud- 
ies were focused on two-dimensional systems or quasi- 
one-dimensional systems with electron injection in the 
lowest subband only. However, in a realistic experimen- 
tal setup the Fermi energy and width of the sample are 
often set at values for which the electron transport is in 
the multi-subband (multi-channel) regime. It is known 
that the SOI will induce interaction between subbands 
as well as the spin statesj 21 ' 22 ' 23 These interactions will 
bring an initially prepared spin state into a quantum- 
coherent spin-mixed state and can give rise to compli- 
cated spin-dependent electron transport phenomena. To 
capture all the physics of the problem an exact multi- 
subband treatment is hence necessary. 

It has been theoretically shown that the vertex cor- 
rections will cancel the transverse spin current in an in- 
finite systemi 12 i 13 i 14 i 15 i 16 However, the disappearance of 
a transverse spin current does not exclude the existence 
of spatial spin density modulations arising as a coher- 
ent interference effect. Reynoso, Usaj and Balseiro have 
recently studied the spin Hall effect as a coherent phe- 
nomena in a finite system by calculating the expectation 
value of the spin operators* 1 ^ In this paper we will calcu- 
late the wavefunction and its individual spin components 
to study the intrinsic spin-Hall effect and Zitterbewegung 
in a multi-channel electron waveguide in the absence of 
an external magnetic field. We will show that the Zitter- 



bewegung and the intrinsic spin Hall effect is essentially 
the same kind of phenomena. We will further show the 
existence of a reversal of the direction of spatially accu- 
mulated spins for electrons injected from individual chan- 
nels in the waveguide, which can not be explained by a 
simple heuristic force operator consideration previously 
adopted in the literature. We will also show that regular 
stripe patterns of spin accumulations are formed inside 
the multi-channel waveguide in addition to the usual edge 
spin accumulations found in the spin Hall effect. These 
patterns are shown to arise from coherent states in the 
waveguide and the number of stripes in the patterns are 
found to be related to the number of open channels in 
the leads. 

To study the problem we consider a two-dimensional 
electron gas (2DEG) formed in the x-y plane of a semi- 
conductor heterostructure. The 2DEG is restricted to a 
waveguide of width w by a hard- wall confinement poten- 
tial U c (y), i.e., U c (y) — for y E [— to/2, to/2] and oo 
otherwise, in the transverse y-direction. The structure 
is subjected to an electrical field along the unit vector £ 
giving rise to a Rashba SOI of strength a. For crystal 
structures without inversion symmetry, there will in ad- 
dition be a Dresselhaus SOI contribution, characterized 
by the strength 0. Such inversion asymmetry is found in, 
e.g., zinc-blende crystals of III-V materials. The waveg- 
uide is connected with two perfect leads with vanishing 
SOI. The electrons are injected from one of the two leads 
into the region with a finite SOI and transport coherently 
in the longitudinal x-direction. The model Hamiltonian 
describing such a system under the effective mass approx- 
imation has the form, 



H = CTr 



2m* 



U c (y) 



1 

2h 



a(r)er • (p X £) 



+ p(r) (a x p x - o-yPy) + B.c. 



(1) 



where p is the momentum operator, m* the electron ef- 
fective mass, which is taken to be to* = 0.042 m e cor- 
responding to an InGaAs quantum-well system, er = 
(a x , o~y , o z ) the vector of the Pauli matrices, and <r the 
unit matrix. The inclusion of the SOI breaks the reflec- 
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tion invariance of the system in the transverse direction, 
i.e., [R y , H] ^ 0. However, for a symmetric confining po- 
tential the system is invariant under operation of & y R y , 
i.e., [a y Ry, H] — 0. This symmetry has an important im- 
plication: for a spin-unpolarized electron injection, the 
spin-up electron probability distribution in the system is 
exactly the mirror image of the spin-down electron prob- 
ability distribution with respect to the transverse reflec- 
tion R y . Thus, if an accumulation of spin occurs on one 
side of the waveguide, the same amount of accumulation 
of the opposite spin will occur on the other side of the 
waveguide. Since an unpolarized electron injection beam 
is an incoherent mixture of two spin injections, this sym- 
metry asserts that all information about the probability 
distribution, regardless of the injection condition, is con- 
tained in the probability distribution of a polarized injec- 
tion beam. One needs only to consider the calculations 
for one spin injection, since the results from the other 
spin injection follows from the symmetry requirement. 

The (spin-unpolarized) probability distribution, 
p J (x,y), for electrons injected from a lead in a spin-7 
state can be obtained from 



(a) 



(b) 



(c) 
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(2) 

where 9^(x,y) is the wave function of an electron in- 
jected from the lead in the qth subband with spin 7 = f 

or I, |x(7=T)> - (M) T and |x(7=i)> = (0,1) T are 
the two spin states with the spin quantization axis along 

1/2 

the electric field direction, k q = [2m * (Ep — E q )/h 2 ~\ 
with Ep being the Fermi energy and E q the qth. subband 
energy in the injection lead, and p c7 ' y (x,y) is the spin-cr 
probability distribution. In the present work, the wave 
function \I/ 7 has been calculated by employing an ex- 
act, spin-dependent, multi-mode scattering matrix tech- 
nique (see Refs. IHiiii for the detail of the theoretical 
method). 

Consider now the case when only the Rashba SOI 
is present, i.e., /3 = and a ^ 0, in Eq. (1). Since 
[tJi,H] ^ for all i, an electron initially prepared in a 
pure spin state in the lead will evolve into a quantum- 
coherently spin- mixed state as it travels in the waveguide 
with a finite SOI. This coherent spin evolution will give 
rise to a spin density modulation along the waveguide. 
To visualize these properties, we calculate the electron 
probability distributions in the waveguide. We consider 
the common situation with the electric field, set along 
the z-direction. In this case, the Rashba SOI term can 
be written as (a/2h)(a x p y — a y p x ) + H.c. 

We first assume the width of the waveguide w = 100 
nm and the Fermi energy Ep = 2 meV, corresponding to 
the case where only one channel is open in the leads. The 
SOI strength is taken to be a = 3 x 10~ n eVm in the 
region with the full strength of the Rashba SOI, defined 
in x e [0, 800] nm, and decreased adiabatically down to 
zero in the transition regions of the entrance and exit, 
defined in x e [-100,0] and x e [800,900] nm, respec- 
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FIG. 1: (Color online) Normalized probability distribution 
for electrons injected in a spin-up state (a)-(b) and for a spin- 
unpolarized electron injection (d)-(f). (a) is for probability 
distribution, p\; (b) for its spin-down projection, pj T ; (c) for 
its spin-up projection, pj T ; (d) for the probability distribution 
of spin-down electrons in the spin-unpolarized injection, pj T + 
p\^; (e) for the probability distribution of spin-up electrons 
in the spin-unpolarized injection, p] T + pj^; and (f) for the 
spin polarization distribution in the spin-unpolarized electron 

T T T I IT II 

injection, p llPO i ~ p[ + p[ ~ Pi Pi ■ in au the figures, 
the channel width is set at w — 100 nm and the Rashba 
SOI strength is a = 3 x 10 -11 eVm for x G [0,800] nm and 
is decreased adiabatically down to zero in x G [—100,0] U 
[800, 900] nm. The Fermi energy is set at E F = 2 meV. The 
upper scale bar shows the measure for plots (a)-(e), while the 
lower scale bar shows the measure for plot (f). 



tively. Figure [Ha) shows the calculated total (charge) 
probability distribution p^ = p^T _|_ pit f or electrons in- 
jected in a pure spin-up state from the lower lead. It 
is seen that the electron wave function exhibits a trans- 
versely oscillating behavior along the waveguide. This is 
the Zitterbewegung arising from interference between the 
two spin components of the electron state. By consider- 
ing the spin-up [Fig QJb)] and the spin-down [Fig QJc)] 
projections separately, spin-wave patterns, in the form of 
spatially separated localized islands of high spin distribu- 
tion density, are found. The two probability distribution 
components, p'^ and p^, for the spin-down injection can 
be obtained from mirror reflection of the results shown in 
Figs.[IJa)-[Hc) as required by [a y R y , H] — 0. It is evident 
that the spin-up and spin-down electron wave-function 
components are localized in different sides of the waveg- 
uide. This accumulation of different spins along opposite 
edges of the waveguide is a fundamental property of the 
spin Hall effect observed recently. The spin projected 
electron probability distributions for a spin-unpolarized 
injection are shown in Figs, [ljd) andQJe). Due to the 
symmetry of the system under the operation of <r y R y , 
the total distribution will be symmetrical in the trans- 
verse direction, and the spin-up (spin-down) probability 
distribution, Fig. [TJd) [Fig. [TJe)], is the mirror image of 
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the corresponding spin-down (spin-up) probability dis- 
tribution. Hence no spin-polarization will be detected in 
the measured two-terminal conductance. This was rig- 
orously shown, in Refs. [26| and [27], to be always hold, 
independent of the details of the conductor, when only 
one channel is open for conduction in the leads. Accu- 
mulations of electron spins with opposite polarizations 
at the opposite edges of a wide conductor have recently 
been experimentally observed&i^. What is measured 
in such experiments is the net spin density distribution, 
Ppoi = p + p — p — p j in the sample for a spin- 
unpolarized electron injection. Figure [TJf ) displays the 
results of calculations for p vo \ in the waveguide and shows 
clearly the signature of the spin Hall effect. It is seen that 
the spin Hall effect and the Zitterbewegung is essentially 
the same kind of phenomena, since one follows from the 
other, depending on if only one or both spin states are 
considered in the injection source and whether the spin 
accumulation or the charge accumulation is measured. In 
other words, the results of Figs. [TJa), [TJd), [TJe) andQJf) 
can be constructed from Figs. QJb) and QJc) and their 
mirror images. 

When the Fermi energy is increased, more channels 
are open for conduction in the leads. Figure [2] shows 
the probability distributions for the Fermi energy set at 
Ep = 12 meV, for which three spin-degenerate channels 
are open for conduction in the leads. In Figs. [2ja)-[2jf), 
the spin probability distributions for spin-up electrons in- 
jected from the lower lead in the three individual channels 
are shown. The probability distributions for the spin-up 
injection in the lowest channel, p\^ and p\\ shown in 
Figs, d^a) and [U^b) exhibit similar Zitterbewegung and 
spin-Hall patterns as was found in the single-channel sys- 
tem. However, here a reversal of the spin-polarization 
direction is found: the spin-up (spin-down) probabil- 
ity density is localized along the left (right) edge of the 
waveguide. This is also true for the probability distribu- 
tions for the spin-up injection in the second lowest chan- 
nel, p| and p'^, as shown in Figs. [He) and^d). Only 
the probability distributions for the spin-up injection in 
the highest subband, pl^ and shown in Figs.^e) and 
(f), exhibit the same spin parity as in the single-channel 
opening case. We note that this phenomenon can not be 
explained by the force operator, derived from the Hamil- 
tonian, as in, e.g., RefllOl. and implies that such a heuris- 
tic approach fails to explain the full quantum mechanical 
calculations. Indeed, for equilibrium systems considered 
here, no physical force is present and the observed Zit- 
terbewegung and spin-Hall effect are non-dynamical but 
merely coherent phenomena. 

In light of these findings, one may ask if the Zitterbe- 
wegung and spin-Hall effect can still be observed in the 
multi-channel system when the contributions from all the 
individual channel injections are taken into account. To 
answer this question we show the total probability distri- 
bution, p\ for a spin-up electron injection in the waveg- 
uide in Fig. [2Kg) and the spin polarization distribution, 
Ppoi, for a spin-unpolarized electron injection in Fig.^h). 
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FIG. 2: (Color online) Normalized electron probability dis- 
tribution and spin-projected electron probability distribution 
for spin-polarized and spin-unpolarized injections in the same 
waveguide as in Fig. [1] at the Fermi energy Ef = 12 meV 
(the case with three opening channels in the leads), (a) is 
for p} T ; (b) for p} T ; (c) for (d) for p^ ■ (e) for pJT. (f) 
for p)jj , (g) is the the normalized total probability distribu- 
tion for spin-up electron injection, p T = X^=i(Pg T + P<j T )i (h) 
is the normalized total spin distribution for spin-unpolarized 
electron injection, p poi = Y%=i(pV + Pg l ~ P? T _ Pq l )\ and 
(i) the same as (h) but for the distribution of the spin polar- 
ization along the transverse y direction. The upper scale bar 
shows the measure for plots (a)-(g) and (i), while the lower 
scale bar shows the measure for plot (h). 



It is seen in Fig. [2tg) that Zitterbewegung oscillations are 
still found along the waveguide, albeit more complicated 
as compared to the single-channel opening case. The 
spin polarization distribution, p po i, shown in Fig. [2fh) 
exhibits a regular spin density pattern. Strong spin po- 
larization with opposite sign is found at the two edges of 
the waveguide, as observed in recent experiments. How- 
ever, spin accumulation, with a pattern of alternative 
spin-polarization stripes, also occurs inside the waveg- 
uide. The number of stripes (including the edge spin po- 
larization stripe) for each spin-polarization direction is 
equal to the number of the open channels in the waveg- 
uide. Experimental observation of the interesting inter- 
nal spin-polarization structure is certainly challenging; it 
requires a spin detection setup with a high spatial reso- 
lution (about 50 nm or better for the system discussed in 
this work). However, it should be noted that no spin Hall 
effect could be observed if the spin polarization along a 
direction perpendicular to the electric field z direction is 
measured in the system. Although the spin polarization 
along the transverse y direction shows a regular stripe 
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FIG. 3: (Color online) Normalized spin probability distribu- 
tion, ppol, for spin-unpolarized electron injection at the Fermi 
energy Ef = 2 meV in a waveguide with (a) w — 370 nm and 
a = 1 x lCT 11 eVm, (b) w = 370 nm and a = 3 x 10 -11 eVm, 
(c) w = 1000 nm and a = 0.1 x 10 -11 eVm, and (d) w = 1000 
nm and a = 3 x 10~ n . The other waveguide parameters 
assumed are the same as in Figs. [T]and[2] 



pattern of accumulation, the same spin polarity is seen 
in all the stripes, as is shown in Fig. [Hi). This can be 
understood as a result of polarization (or magnetization) 
by the effective transverse magnetic field arising from the 
Rashba term. 

For a wide conductor, a large number of channels are 
typically open for conduction and the inter-subband mix- 
ing induced by the SOI can produce a complicated pat- 
tern in the electron probability distribution. In Fig. [3] 
we show the calculated spin polarization p po i for a spin- 
unpolarized electron injection at the Fermi energy Ep = 
2 meV in a waveguide of width w = 370 nm (with five 
open channels) and a waveguide of width w = 1 /im (with 
14 open channels) at weak and strong SOI strengths. The 
distributions are symmetric with respect to <J y R y . For 
the weak SOI strength, the same structural patterns of 
alternative spin-polarization stripes [Figs [3]Ja) and HD^c)] 
as observed in the narrower waveguide are found. For the 



strong SOI strength, the spin-polarization distribution 
patterns [Figs. [3{b) and[3fd)] show complex structures. 
Although spatial spin accumulations are observable even 
under these circumstances, regular spin accumulations at 
the edges of the waveguides and internal stripe structures 
are destroyed by the strong SOI induced inter-channel 
scattering. This result indicates that the regular spin 
Hall effect can only be clearly observed in a range of SOI 
strengths in a wide, multi-channel electron waveguide. 

Since the Rashba and the Dresselhaus term in the 
Hamiltonian are related through a unitary transforma- 
tion, with the substitutions of a x — > a y , a y — > a x and 
z — > — z, the fundamental results presented above hold 
for the case of a = and j3 ^ 0, i.e., when only the Dres- 
selhaus SOI is present. For a = ±(3 the spin-dependent 
part of the Hamiltonian takes the form 

(X 

H± = - {p y ± p x ) er • (e x =F e y ) . (3) 

The SOI induced effective magnetic field direction is in 
this case independent of the momentum. This means 
that a well-defined spin quantization axis can be found 
throughout the SOI region and spin is therefore a con- 
served quantity. In Ref. 28, it was shown that in III-V 
systems the Rashba and Dresselhaus SOIs can be of com- 
parable strengths. This would imply that no Zitterbewe- 
gungi& nor spin-Hall effect could occur in this situation. 

In conclusion, we have studied the total and spin- 
resolved electron probability distributions in a ballistic 
waveguide with SOI. For a spin-polarized electron injec- 
tion, a Zitterbewegung pattern is observed in the proba- 
bility distribution of the electrons in the waveguide. For 
a spin-unpolarized injection, the spin-resolved electron 
probability in the waveguide shows spin accumulations. 
The intrinsic spin Hall effect, namely accumulations of 
opposite spins at the lateral edges of the waveguide, can 
be observed in the SOI region. We have also predicted 
the existence of a regular stripe pattern of spin accu- 
mulations inside the waveguide. Finally we have shown 
that the Zitterbewegung and the spin Hall effect found in 
this work stem from the same mechanism and are formed 
from the coherent states of electrons in the waveguide. 
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